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Abstract: We review the Sakai-Sugimoto model of holographic QCD at zero temperature 
and finite chemical potential, comparing the results to those expected at large- N c QCD, and 
those in a closely related holographic model. We find that as the baryon chemical potential 
is increased above a critical value, there is a phase transition to a nuclear matter phase, the 
details of which depend on the model. We argue that the nuclear matter phase is necessarily 
inhomogeneous to arbitrarily high density, which suggests an explanation of the "chiral 
density wave" instability of the quark Fermi surface in large- N c QCD. Some details of the 
instanton distribution in the holographic dual are reminiscent of a Fermi surface. This short 
manuscript summarizes a talk given by M.R. at "Theory Canada 4" conference, and is based 
largely (but not entirely) on the results of [1]. 
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1. Introduction 

The phase diagram of QCD as a function of temperature and baryon chemical potential displays 
a rich variety of phases. A better understanding of the details of the phase diagram at intermediate 
chemical potentials would have valuable applications, for example in understanding the physics of 
neutron-star interiors. Currently our knowledge of that regime of the phase diagram is based exclu- 
sively on extrapolations and semi-empirical toy models. However, with the advent of gauge-gravity 
dualities, we have a new tool for studying the properties of certain strongly coupled gauge theories. 
It is interesting to study these QCD-like theories in regimes for which neither analytic or numerical 
studies are currently possible in real QCD. In the present paper, we will focus then on equilibrium 
properties of those holographic models, at zero temperature and finite baryon chemical potential. 

Additionally, the study of the phase diagram of QCD involves concepts common in many con- 
densed matter systems, such as Fermi surfaces and their instabilities, superconductivity, and inhomo- 
geneous phases. Therefore, the holographic study of the phase diagram is a natural first step towards 
a better understanding of holographic duals of these interesting phenomena in the context of strongly 
interacting condensed matter systems. 

In order to get behavior similar to real QCD in the regime discussed here, it seems essential to 
study a theory with baryon charge carried exclusively by fermionic fields, otherwise the physics at 
finite chemical potential is likely to involve Bose condensation rather than the formation of a Fermi 
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surface. Such a model was constructed a few years ago by Sakai and Sugimoto [2], and will be focused 
on here. In addition to reviewing the results in [1], we also provide a point of comparison by presenting 
results in parallel for a similar model, defined below, which has flavoured boson-fermion degeneracy 
and is therefore expected to have different qualitative behavior. We present the qualitative results of 
this comparison below. 

2. Overview of Models 

2.1. The Adjoint Sector 

The basic setup for our models begins with the low-energy decoupling limit of N c D4-branes 
wrapped on a circle of length 2irR, with anti-periodic boundary conditions for the fermions. Apart 
from N c , this theory has a single dimensionless parameter A, the four-dimensional 't Hooft coupling at 
the Kaluza-Klein scale. As pointed out by Witten [3], at small A the dynamics should be that of pure 
Yang-Mills theory. Unfortunately, in this limit, the dual gravity background is highly curved, and so we 
are not in a position to study it perturbatively. For large A, on the other hand, the gravity background 
is weakly curved, and so studying the theory is possible by using classical gravity calculations. This 
is a general phenomenon in the gauge-gravity duality; the gravity dual necessarily deviates from an 
asymptotically free field theory in the ultraviolet, and can be thought of as strong (bare) coupling 
expansion. 

In the large-7V c limit the gravity dual is the near-horizon geometry of the D4-branes. The metric 
is [3, 4] 

3 3 

ds 2 = (-^) 2 (v^d*" + f{U)dxl) + 2 (j^ + U 2 dnj) , (2.1a) 

W) = l- (jf) . (2.1b) 

where dQ.4 is the volume-element of a unit 4-sphere. The 24 direction, corresponding to the Kaluza- 
Klein direction in the field theory, is taken to be asymptotically periodic, with coordinate periodicity 
2ttR. It combines with the radial coordinate u to form a cigar-like geometry, with the topologically of 
a disk. 

The study of the pure Yang-Mills theory holographically is similar to AdS/CFT; for example, one 
has the following qualitative features: 

• Harmonics of closed string fields represent the glueball spectrum, while those of open strings, 
introduced below, are related to the meson spectrum. 

• There is a deconfinement phase transition, related to the fact that there are two Euclidean ge- 
ometries contributing to the finite temperature partition function. In fact, those two geometries 
are related by a "flop" - an exchange of the compact direction x 4 with the Euclidean time. 



2.2. Adding Flavours 

Now that we have defined the adjoint sector of the theory, we would like to add Nf flavours of 
fundamental quarks. By keeping Nf fixed in the large-A^ limit, we work in the so-called "quenched 
approximation", whereby the addition of matter (in the form of probe branes) does not modify the 
background geometry of the gravity dual. 

The Sakai-Sugimoto construction is motivated by the observation that the light open string modes 
which live at the (3 + 1) -dimensional intersection of a particular configuration of D4- and D8-branes 
give rise to chiral fermion fields on the intersection without accompanying bosons. Thus, to the original 
D4-branes, which we can take to lie in the 01234 directions, one adds a stack of Nf D8-branes forming 
a line in the cigar spanned by (U, £4). The configuration is represented by the following table: 



©2008 NRC Canada 



M. Rozali, L. Brits 





a; x 1 x 2 x 3 (x 4 ) 


u n 4 


D4: 
D8: 


— — — — • 


• • 



Fig. 1. Brane configuration of the Sakai-Sugimoto model 

The specific embedding of the D8-branes in the bulk depends on the asymptotic separation of the 
stacks (and also on any distribution of matter on the branes), but we will focus exclusively on the 
case where the two asymptotic parts of the D8-brane stack sit at opposite sides of the D8 circle, in 
which case each side simply extends to the tip of the cigar along a line of constant x\. Unlike the 
adjoint fields, the fundamental fields, which arise from 4-8 strings, do not start out supersymmetric, 
and include massless fermions along with heavy bosons. 

Let us compare this case to a different theory, where supersymmetry is not broken in the flavour 
sector. To describe the embedding of the probe brane in that case, let us first introduce new coordinates 
which are more adapted to supersymmetry. To this end, introduce a coordinate p which makes second 
term in Eqn. (2.1a) manifestly conformally Euclidean: We choose a reparameterization u{p) that solves 

* = = P \ (2.2) 

9pp 

i.e., 

2 2 

u \ / i i 



(«) 



which leads to the form of the metric in [4]. The background metric (Eqn. 2.1a) can now be written 
entirely in terms of p, as 

3 3 1 

ds 2 = [£) * (j + {p) Vltu dardar + J^dxl) + R!U p i P) ' 2 (dp 2 + P 2 dn 2 ) , (2.4) 
where we have defined the following functions for convenience: 

/ ± (ri = ( 1± H) = (i^). <«> 

Finally, let us refactor the conformally Euclidean R 5 into R 1 x R 4 , letting p 2 = A 2 + r 2 , so that 

dp 2 + p 2 <m\ = dr 2 + (d\ 2 + X 2 dQl) . (2.6) 

In these coordinates, the embedding of the probe branes in the geometry is described by the following 
table: 

x° x 1 x 2 x 3 (x 4 ) r A Q 3 



D4: 
D8: 

Fig. 2. Brane configuration with boson-fermion degeneracy 

Compared to the Sakai-Sugimoto model, the adjoint sector is unaltered, while the fundamental 
sector now preserves N = 1 supersymmetry (before compactification). Both fundamental scalars and 
fermions now acquire mass of order Mkk- The theory is pure YM at sufficiently low energies, and 
becomes approximately supersymmetric at the KK scale. This is a natural theory to compare to the 
Sakai-Sugimoto model: on the one hand the gravity duals are fairly similar to each other, but on the 
other hand the field theories, especially at finite chemical potential, are expected to be radically differ- 
ent due to the difference in both the massless and the light (QCD scale) matter content. 
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2.3. Baryons 

Configurations with non-zero baryon charge density require sources for Aq, the gauge field on the 
D8-branes. The basic source for A is the endpoint of a fundamental string. In order to have some net 
charge, we need a source for fundamental strings in the bulk. In our background we can get a density of 
charge on the D8-brane by having a density of D4-branes wrapped on S 4 in the bulk, with iV c strings 
stretching between each D4-brane and the D8-brane. In the case where we have Nf > 1 D8-branes, 
our baryon branes can dissolve in the D8-branes and show up as instantons. The question of which of 
these two pictures is more appropriate is a dynamical one. 

In order to work at finite chemical potential [i for baryon number, we need to turn on non-normalizable 
mode for the operator dual to the baryon number. We require then that the value of the gauge potential 
A approaches the constant [i in both asymptotic regions of the D8-brane. 

3. One Flavour Physics 

In this section we study the physics of both models at finite chemical potential, in the simpler case 
of a single quark flavour. Since we have only a single D8-brane in the bulk, we can use the Abelian 
Born-Infeld action for our analysis. Baryons are represented by pointlike charges on the D8-branes 
that source the electrostatic potential; configurations with smooth baryon density correspond to having 
some density of these pointlike instantons on the D8-brane. 

The dynamics of the probe D8-branes are governed by the action: 



Sdbi = -Ms J d 9 <re 4 ' v / - dct(g ab + 2na'F ab ) 

Here g a b is the induced metric on the brane, and F a b is the field strength on the brane's worldvolume. 
There is also a Wess-Zumino term, which will be irrelevant for the discussion. 



3.1. The Sakai-Sugimoto Model 

Let us discuss the Sakai-Sugimoto model first. To simplify the action, we can work in the "static 
gauge", by identifying the worldvolume and spacetime coordinates in the sphere and the field theory 
directions, as well as in the radial direction U. Additionally we choose to concentrate on the case where 
X4 is constant. We are interested only in time-independent configurations which are homogeneous 
and isotropic in the spatial directions of the field theory (labeled by indices i, j, k). For details of the 
derivation see [1]. 

For any value of the chemical potential, which is the asymptotic value of the gauge potential A n , 
the constant configuration A (X) = M is always a solution of the equations of motion. For this solution 
there is no electric flux, and by Gauss law we also have zero baryon density. Whether or not other solu- 
tions exist, with some non-zero baryon density, depends on the value of \i. This leads to the existence 
of distinct phases as function of [i and phase transitions between them. 

Once we have some non-zero baryon density, in our setup the embedding of the branes does not 
change. However, to find the configuration of the gauge field Aq it is important to include the energy 
contributions from the electric field on the branes. Starting with the DBI action above, we can derive 
the 3+1 dimensional energy density via a Legendre transformation. We find 

c Ms 8 2 | f U§ 

tfiux = 7r -ft? / du — . . (3.7) 

1 9s 3 4 J v //(l-2W/((7)(a l ;Ao) 2 ) 

As a first approximation, we make the simplifying assumption that all the pointlike instantons sit 
at U = Uq. In this simple approximation, the energy from the electric flux E — S ^ A ) (from both 
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halves of the D8-brane) is 



^8 16^ 2 

9s 3 



U 

Tr 2 RjU£h(e), (3.8) 



where we have subtracted off the E = contribution from the brane tension, and defined e = E/Uq 
and 



/>oo 
dx(\/ 



h(e) — / dx{y x 5 + e 2 — x 2 ) 



1 



Meanwhile, the energy from the charges in the electrostatic potential and the masses of the pointlike 
instantons combine to give a term 

^charge = — (M ~ Mc)«B • 

where [i c is a positive constant. 

Defining ft = ^fff and using the relation between ub and E (which is the DBI version of Gauss 
law, relating the charge to the electric flux), the total energy may be written as 

M8 16 2 D f rr | f u ,^ 1, 



£ = f s ^ 2 Rm\h{e)--(n-i)ey 

From this, we find that there is a non-trivial solution provided /j, > fj, c , the total energy is then mini- 
mized when 

±(H-l) = h'(e). 

This can be inverted to determine the relationship between n b and [i above the transition, we find that 
the chemical potential exhibits the following asymptotic behavior: 

ub oc fi — [i c small yU — fi c , (3.9) 

5 

n B oc ^2 large /U — /i c . (3.10) 

In summary, for a given value of the chemical potential greater than the critical value, we find 
some preferred distribution of charges on the D8-brane. The total baryon density for a given value of 
li may be read off from the asymptotic value of the electric flux, and the result increases smoothly 

5 

from above the critical chemical potential, approaching an asymptotic behavior tib oc pa . With a 
more careful analysis, allowing for dynamical charge distribution, it was found [1] that the qualitative 
behavior of ns{p) is the same as in the simplified model of the previous section, though the numerical 
coefficients come out different. 

Additionally, the behavior of the energy £{p) can be found 

£ oc {p — pc) 4 small \i — yU c , (3.11) 

7 

£ oc fi? large p — /j, c . (3.12) 

We see that the energy vanishes at the critical point, signaling a second order phase transition. For the 
single flavour case, even allowing for dynamical charge distribution, the transition to nuclear matter 
is continuous, unlike QCD, but it may be expected that the single flavour case is different due to 
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the absence of pions which usually play a crucial role in interactions between nucleons. Indeed, the 
details found in [1] are different for the multi-flavour case, the transition is first order, and the critical 
exponents at large p are also different 



n B oc ji 3 (3.13) 



£ oc p 4 (3.14) 

Note that these are the powers appropriate for free fermions, and we would like to understand this point 
better. 

3.2. Supersymmetric Matter Content 

Let us now analogously develop the second model, in which there is a Bose-Fermi degeneracy in 
the fundamental matter. For the worldvolume coordinates of the D8-brane we choose the "static gauge" 
except for r(A). Thus, the worldvolume coordinates are {x°...x 4 , A, The induced metric is on the 
D8-brane is 

dsi = (-£) ' (f+ip^dx^dx" + j^dx?) + [(1 + d x r 2 )d\ 2 + A 2 ^ 2 ] . (3.15) 

Let us allow one component of the gauge potential on the D8 world-volume to be non-zero, in response 
to a source of U(l) charge that is concentrated at a point U = IV Let A = Ao(X) so that only 
F\o = d\Ao is non-zero. We need to evaluate the following determinant: 

det(g + F) = _/ + (p)V_(p) 2 ^ \-d x A (\) 2 + f + (p)Hl + dxr(\) 2 )] (dQ 3 ) 2 . (3.16) 

The DBI action governing the D8 world-volume is therefore given by 

2tt 2 (j,&V 5 



S 8 = - 



J rfA/+(p)/_(p)A 3 [-d x A (\) 2 + /+(p)t(l + d\r(X) 2 ) 



(3.17) 



This system does not admit a trivial (constant) solution in the coordinates that we have chosen, and 
so we will require a careful analysis of the boundary conditions due to the addition of baryons in order 
to determine the most favourable embedding. Let the number density ps(A) = ubS(X) represent a 
localized source of baryons. We minimally couple baryons to the gauge field via addition of the term 

S charge = N c jd\A (\) PB (\), (3.18) 

p CO 

= N c n B p-N c n B / d\d x A a (X). (3.19) 
Jo 

The kinetic contribution of the baryons is [1] 

S m ass = ~^-^7r 2 Rln B f d\U(X)5(X) = -n B N c M B . (3.20) 
9s «J J 

To find the equations of motion, we note that the variable A is effectively cyclic, except at the 
boundary, so that by varying A, we find again that the flux 

P SS 2n 2 p 8 / + (p)/_(p)A 3 

E = = d x A (A), (3.21) 

A 9s ^-d x A (\)* + f + (p)t(l + d x r(\) 2 ) 
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is conserved (more precisely, E — N c ub)- Thus we have the equation of motion 

d x Ao(X)=E 



/ + (p)i(l + A r(A)2) (322) 



Roughly this corresponds to the following asymptotics as A — > oo: 

where E 1 = HbN c [5]. To simplify calculations, we perform a Legendre transform with respect to 
dxA (X): Define 

S* S =S 8 - Jd\Ed x A (E), 

= ^B. ( dX f + (p) I y/l + dxrWy/MpFf-WXe + E*. (3.24) 

where we have defined 

E=*E = ?- = !!pL. (3.25) 
The equations of motion for r(A) follow from the action 

q* rr* , Q 

- S 8 * - — ^7r 2 i?3 ns /" d \ U{\)6{\) (3.26) 

To obtain the boundary condition on d\r(X = 0), we formally extend A to be negative and extremise 
the "doubled" system. Then 

d ( SS* \ SS* n 

= 0. (3.27) 



d\ \8d x r(X) J Sr(X) 
We may integrate from — e to +e and find 

2 #fr - °W + ~^ 2r >b [ +e dX ^(X)S(X) = 0, (3.28) 
ddxr(e) g s 3 y_ e dr 

and, taking the limit e — > 0, we obtain the required boundary condition: 

2 ^S) = -ab^Se* (3 - 29) 

9 2 

where fcg = x M8 (2ttR). We may use this to solve for 9at(0) in terms of r(0) and ng. Here it is 
important to note that ^ is strictly negative, so that the solution consistent with Eqn. (3.29) must be 
selected. With the boundary condition determined, we may ignore the baryon and solve the homoge- 
neous equations of motion for r(X) numerically. After specifying r(oo) as the bare quark mass, we can 
then solve for r(0) in terms of r(oo). Each embedding is therefore uniquely specified by the additional 
parameters [i and E, or, equivalently, p and ns- 
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The relationship between ub and p is fixed by requiring that the total energy of the configuration 
be minimized. Since the configuration is static, the energy is equal to minus the Lagrangian evaluated 
over a solution to the equations of motion, i.e., 

£ = — jy (5*8 + S c harge + S mass ) , (3.30) 

where V4 — Jdx° J d 3 x. For numerical purposes, we work in string units, noting that the number of 
integrations eliminate all powers of £ s . The first two terms partly combine to form 

£ f h* = hj d\f + { P )i^\ + d x r{\f^f + {pyf_( P y\s + E^ (3.3i) 

while the remainder is given by —p,N c riB and £ maS s = ^U m i n riB — ubN c Mb. All contributions 
together give the energy in terms of the convenient variable E, which is then to be minimized over the 
family of solutions r(A) parameterized by this variable. 

f = - (p - Mb) E + f°°d\ U (p) 1 y/l + dxr(\)*Jf + (p)*f-(p)*\* + E2. (3.32) 
Jo 

For a fixed value of p, we may seek the value of E which solves 

= dl { Mb ^ + P A/+(p)lv/i + d ^ X ) 2 \/f+(p) 2 f-(p) 2X6 + E2 ) ■ ( 3 - 33 ) 

As it stands, Eqn. gives the chemical potential in terms of baryon number, without modification. All 
quantities may be numerically evaluated, and the behaviour of p(ns) for small values of ns is shown 

2 

in Fig. (3.2), along with a best fit to p, — )jl c ~ n B , with p c w 0.32 Mb(E = 0). 

1.4 

1.2 
1.0 
0.8 
0.6 
0.4 
0.2 

0.1 0.2 0.3 0.4 

Fig. 3. Power-law dependence of chemical potential on baryon number with best fit indicated by solid line. Shown 

Too = |po- 

The calculation presented so far has assumed a non-zero flux E. There is always another solution, 
with E = 0, where p may take on any value, corresponding to a constant gauge potential Aq(X) = p 
in the absence of sources. For values p < p c , this solution is energetically favourable, after which 
baryons begin to accumulate. This situation has been analyzed using the best-fit data obtained above, 
which indicates a smooth transition of second or higher order into the nuclear matter phase. 

It is interesting to compare the behavior here to that found in the Sakai-Sugimoto model. We have 
found that the details of the transition to nuclear matter are identical; it is smooth in both cases, with 
same critical exponent near the transition point. This is somewhat surprising given that the massless 
matter content is very different in both cases. Perhaps the Nt = 1 case cannot be expected to be very 
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different from pure YM, and quanlitative differences only arise when massless pions become part of 
the story. 

On the other hand, we do find different critical exponents at large chemical potential. This is ex- 
pected, as in our model there is an approximate fermion-boson degeneracy at high enough energy, 
which is sufficient to cause the model to behave differently in the regime of high chemical potential. 

Most interestingly, we find that the binding energy per nucleon is of the order of magnitude it is 
expected to be, the baryon mass itself. In other words the baryons form a deep bound state, unlike the 
case of the Sakai-Sugimoto model, and of real-world QCD. It is interesting to understand further the 
mechanism for forming shallow bound states, and the dependence of the binding energy on specifics 
of the model. 

4. Other Results for the Sakai-Sugimoto Model 

We review here other results obtained in [1], it would be interesting to compare those to other 
models, to probe their dependence on the precise matter content and other details of the theory. 

4.1. Transition to Nuclear Matter 

We have seen that turning on a chemical potential in the gauge theory corresponds to including 
boundary conditions Aq = fi for the two asymptotic regions of the D8-brane. For any [i, one solu- 
tion consistent with these boundary conditions is to have constant Aq everywhere on the brane. This 
represents the vacuum configuration in the field theory. However, beyond a certain critical chemical 
potential, this solution is unstable to the condensation of baryons. This is manifested by the existence 
of another classical solution with non-trivial flux, which dominates the thermodynamics beyond the 
critical value /i c . A first approximation to this critical value is the baryon mass, since it is at this point 
where it becomes energetically favorable to add single baryons to the vacuum. In fact, the critical value 
is somewhat lower, since the baryons have binding energy. Therefore the critical value /i c can be used 
to estimate the binding energy per baryon. 

The binding energy per nucleon is estimated in real QCD to be 

b vol = 16 MeV . (4.34) 

which is notably small compared to the QCD scale. 

In the Sakai-Sugimoto model, for a single flavour the transition is second order, thus involving no 
latent energy. Within the localized source approximation the critical point is precisely at the baryon 
mass /i c = Mb, therefore in that case the binding energy per nucleon vanishes. However, in the case 
with Nf > 1, we can have nonsingular instantons on the Nf coincident D8-branes, and the minimum 
energy configurations for large enough [i are should involve smooth configurations of the non- Abelian 
gauge field carrying an instanton density. This is the mechanism by which the baryons bind to each 
other in the holographic context. 

For many flavours, the Sakai-Sugimoto model displays a first order transition to nuclear matter at 
some calculable critical chemical potential yL c . We find the binding energy per nucleon 

Ebind = M B - ii c i=a ^M K k(c' - c) , 

where c, c' are unknown constants of order one, and the result is insensitive to the value of A for large 
A. Since we also know that this binding energy approaches some constant value in the limit of small 
A, then assuming a smooth behavior at intermediate values of A, we can treat the large A result as a 
prediction for the order of magnitude of the QCD result. Noting that M KK Aq C d for large A, the 
value of the binding energy per nucleon extrapolated to N c = 3 becomes 

Ebind = ^-A QCD (c' - c) w 7 MeV(c' - c) 
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thus we obtain the same order of magnitude as the QCD result. In particular the binding energy is 
anomalously small compared to the QCD scale. It would be interesting to understand the mechanism 
for this, and how sensitive is the result to exact details of the model. 

4.2. Chiral Density Waves 

One significant difference between the large-iV c theory and ordinary QCD is the expected behavior 
at asymptotically large values of the chemical potential [6]. In both cases, we have attractive interac- 
tions between excitations on the Fermi surface that result in an instability, but the nature of the resulting 
condensates is different. Whereas for N c = 3 the instability is a BCS-type instability, believed to lead 
to a color superconductor phase, the dominant instability at large-iV c is toward the formation of "chi- 
ral density waves", inhomogeneous perturbations in the chiral condensate with wave number of order 
twice the chemical potential. This suggests that the ground state for large-7V c QCD at large enough 
chemical potential is inhomogeneous, however the nature of the true ground state remains mysterious. 

In the holographic model of Sakai and Sugimoto, it was found in [1] that there are no allowed 
configurations of the D8-brane gauge field that are spatially homogeneous in the three field theory di- 
rections such that the net energy density and baryon density in the field theory are both finite. Thus, 
any phase with finite baryon density is necessarily spatially inhomogeneous. This has a simple inter- 
pretation: it suggests that at large-7V c , the nucleons retain their individual identities for any value of the 
chemical potential. This suggests that the chiral density wave instability of the quark Fermi surface in 
large- N c QCD simply indicates that the quarks want to bind into nucleons even at asymptotically large 
densities. It is interesting to understand this phenomena further. 

4.3. Holographic Fermi Surface 

It is interesting to observe some characteristics of the charge distribution, once we allow it to 
become dynamical. This is summarized by a density of pointlike instantons, in the single flavour case, 
and by the distribution of the instanton density in the non-Abelian case. In all those cases, the charge 
distribution has a sharp edge in the holographic radial direction, at U = U max which progresses 
further and further towards the UV as the chemical potential is increased. In the field theory picture, 
the radial direction represents an energy scale, so the charge distribution we find in the bulk should be 
related to the spectrum of energies for the condensed baryons. The edge of the distribution may be a 
bulk manifestation of a (quark) Fermi surface. It would be interesting to check this interpretation by 
comparing to other models with flavoured massless bosons, where such Fermi surface is not expected. 
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